Auditory sensitivity provided by 
self-tuned critical oscillations of hair cells 
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We introduce the concept of self-tuned criticality as 
a general mechanism for signal detection in sensory 
systems. In the case of hearing, we argue that active 
amplification of faint sounds is provided by a dynam- 
ical system which is maintained at the threshold of 
an oscillatory instability. This concept can account 
for the exquisite sensitivity of the auditory system 
and its wide dynamic range, as well as its capacity to 
respond selectively to different frequencies. A spe- 
cific model of sound detection by the hair cells of the 
inner ear is discussed. We show that a collection of 
motor proteins within a hair bundle can generate os- 
cillations at a frequency which depends on the elastic 
properties of the bundle. Simple variation of bun- 
dle geometry gives rise to hair cells with character- 
istic frequencies which span the range of audibility. 
Tension-gated transduction channels, which primar- 
ily serve to detect the motion of a hair bundle, also 
tune each cell by admitting ions which regulate the 
motor protein activity. By controlling the bundle's 
propensity to oscillate, this feedback automatically 
maintains the system in the operating regime where 
it is most sensitive to sinusoidal stimuli. The model 
explains how hair cells can detect sounds which carry 
less energy than the background noise. 

Detecting the sounds of the outside world imposes 
stringent demands on the design of the inner ear, where 
the transduction of acoustic stimuli to electrical signals 
takes place The hair cells within the cochlea, which 
act as mechanosensors, must each be responsive to a 
particular frequency component of the auditory input. 
Moreover, these sensors need the utmost sensitivity, since 
the weakest audible sounds impart an energy, per cycle of 
oscillation, which is no greater than that of thermal noise 
At the same time, they must operate over a wide 
range of volumes, responding and adapting to intensities 
which vary by many orders of magnitude. Clearly, some 
form of non-linear amplification is necessary in sound de- 
tection. The familiar resonant gain of a passive elastic 
system is far from sufficient for the required demands, be- 
cause of the heavy viscous damping at microscopic scales 
Instead, the cochlea has developed active amplifica- 
tory processes, whose precise nature remains to be dis- 
covered. 

There is strong evidence that the cochlea contains 
force-generating dynamical systems which are capable of 



executing oscillations of a characteristic frequency j§- 10 1 . 
In general, such a system exhibits a Hopf bifurcation [ 1 1 1 : 
as the value of a control parameter is varied, the behavior 
abruptly changes from a quiescent state to self-sustained 
oscillations. When the system is in the immediate vicin- 
ity of the bifurcation, it can act as a nonlinear amplifier 
for sinusoidal stimuli close to the characteristic frequency. 
That such a phenomenon might occur in hearing was first 
proposed by Gold ||, more than 50 years ago. The idea 
was recently revived by Choe, Magnasco and Hudspeth 
]l2[ , in the context of a specific model of the hair cell. No 
general analysis of the amplification afforded by a Hopf 
bifurcation has been provided, however, and no theory 
has been advanced to explain how proximity to the bi- 
furcation point might be ensured. 

In this paper, we provide both a generic framework 
which describes the known features of acoustic detec- 
tion, and a detailed discussion of the specific elements 
which could be involved in this detection. We first derive 
the general resonance and amplification behavior of a dy- 
namical system operating close to a Hopf bifurcation and 
emphasize that such a system is well-suited to the ear's 
needs. In order for active amplification to work reliably, 
tuning to the bifurcation point is crucial. We introduce 
the concept of a self-tuned Hopf bifurcation which per- 
mits the favorable amplificatory properties of a dynami- 
cal instability to be obtained in a robust way. Self-tuning 
maintains the system in the proximity of the critical point 
and is achieved by an appropriate feedback mechanism 
which couples the output signal to the control parameter 
that triggers the bifurcation. The concept can explain 
several important features of the auditory sensor such as 
the frequency selectivity, high sensitivity and the ability 
to respond to a wide range of amplitudes. It can also 
explain the intrinsic nonlinear nature of sound detection 
@0 and the occurrence of spontaneous sound emission 
by the inner ear Furthermore, self-tuned criticality 

provides a framework for understanding the role of noise 
in the detection mechanism. The amplificatory process, 
which involves a limited number of active elements, intro- 
duces stochastic fluctuations, which adds to those caused 
by Brownian motion. We show that the response to weak 
stimuli can take advantage of this background activity. 

The proposed existence of a self-tuned Hopf bifurcation 
raises questions about the specific mechanisms involved: 
What is the physical basis of the dynamical system? 
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How is the self-tuning realized? It might be expected 
that different organisms have evolved different apparatus 
to implement the same general strategy. In this paper, 
we restrict our specific discussion to the more primitive 
cochleae of non-mammalian vertebrates. We propose a 
model of the hair cell of the inner ear which accords with 
data from a wide variety of physiological experiments. 
The model incorporates a physical mechanism which al- 
lows motor proteins to generate spontaneous oscillations 
p5[ . We find that molecular motors such as dyneins in 
the kinocilium or myosins in the stereocilia are natural 
candidates for the force generators involved in the ampli- 
fication of hair-bundle motion. Tension-gated transduc- 
tion channels in the stereocilia serve primarily to detect 
this motion, but also have a second function: by admit- 
ting ions which regulate the motor protein activity, they 
provide the self-tuning mechanism. 



I. GENERIC ASPECTS 

Amplification and frequency filtering of a Hopf 
bifurcation. We discuss the behavior of a dynamical 
system which is controlled by a parameter C. Above a 
critical value, C > C c , the system is stable; for C < C c it 
oscillates spontaneously. At the critical point (or Hopf bi- 
furcation) C — C c the system shows remarkable response 
and amplification properties which do not depend on the 
physical mechanism at the basis of the bifurcation. These 
generic properties can be described as follows. Since 
we are interested in the response to a periodic stimu- 
lus with frequency v = u/27r, we express the hair-bundle 
deflection x(t) by a Fourier series x(t) = x n e mult with 
complex amplitudes x n = x*_ n . In the vicinity of the 
bifurcation, the mode n = ±1 is dominant and the re- 
sponse to an externally-applied sinusoidal stimulus force 
f(t) = fie tuJt + f-ie~ luJt can be expressed in terms of 
a systematic expansion in x\ . Symmetry arguments (see 
Appendix A) imply that the first nonlinear term is cubic: 



h =A Xl +B 



(1) 



where A(u>, C) and B(u>, C) are two complex functions. 
The bifurcation point is characterized by the fact that A 
vanishes for the critical frequency, A(uj c , C c ) = 0. For 
C < C c and no external force, the system oscillates with 
|a;i| 2 ~ A 2 (C c — C)/C c , where A is a characteristic 
amplitude. For C = C' c the response to a stimulus at the 
critical frequency has amplitude 



\B\ 
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This represents an amplified response at the critical fre- 
quency with a gain 



l/ll 



(3) 



that becomes arbitrarily large for small forces. 

If the stimulus frequency differs from the critical fre- 
quency, the linear term in Eq. (|l|) is non-zero and can be 
expressed to first order as A(cu, C c ) ~ Ai (to — lu c ). The 
dramatic amplification of weak signals, implied by Eq. 
(||), is maintained as long as this term does not exceed 
the cubic term in Eq. (|l]). If the frequency mismatch 
increases such that \lj - ui c \ > l/il 2 / 3 !^ 1 / 3 /^!, the re- 
sponse becomes linear 



l/i 



\(UJ - LU c )Al\ 



(4) 



i.e. the gain is independent of the strength of the stimu- 
lus. 
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FIG. 1. Response to external forces near a Hopf bifurca- 
tion (a) Amplitude |cci| as a function of force at various 
driving frequencies u (O 2 kHz, x 5 kHz, □ 10 kHz, + 13 
kHz), (b) Gain r as a function of frequency v for different 
amplitudes (O 0.01 pN, A 0.05 pN, + 0.1 pN, x 0.5 pN, 
□ 1 pN). While the form of these curves is generic, the numer- 
ical values of force, amplitude and frequency depend on the 
physical nature of the dynamical system. The values given 
here correspond to the specific model model of a hair cell dis- 
cussed in the text, with parameters chosen to give a critical 
frequency of approximately 10 kHz. 

Thus the Hopf resonance acts as a sharply tuned high- 
gain amplifier for weak stimuli, and as a low-gain filter 
for strong stimuli. This generic behavior is illustrated 
in Fig. 1 with data obtained by numerical simulation. 
Laser interferometry measurements of the motion of the 
basilar membrane, when the live cochlea is stimulated 
by pure tones, display strikingly similar features 
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In particular, the peak response as a function of force 
amplitude has been demonstrated to obey a power law 
M ~ |/i|°' 4±0 ' 2 @- This strongly suggests that the 
membrane is being driven by a dynamical system which 
is poised at a Hopf bifurcation. 

Self-tuned critical oscillations. How does the sys- 
tem come to be so precisely balanced at the critical 
point ?Q The control parameter must be tuned to C ~ C c , 
otherwise the nonlinear amplification is lost. Moreover, 
the value of C c differs for hair cells with different char- 
acteristic frequency. We propose a feedback mechanism 
which allows the dynamical system to operate automat- 
ically close to the bifurcation point, whatever its char- 
acteristic frequency. Without loss of generality, we as- 
sume that the control parameter decreases as long as the 
system does not oscillate. After some time, critical con- 
ditions are reached and spontaneous oscillations ensue. 
The onset of oscillations triggers an increase of the con- 
trol parameter which tends to restore stability. Hence 
the system converges to an operating point close to the 
bifurcation point. To illustrate this general idea, we con- 
sider the following simple feedback which changes C in 
response to deflections a;: 



]_dC _ 1 (x*_ 
C dt ~ t U 2 



(5) 



where S is a typical amplitude. If no external force is ap- 
plied, this feedback, after a relaxation time r, tunes the 
control parameter to a value Cs for which spontaneous 
oscillations with \xi \ ~ S occur. If 6 is small compared to 
the characteristic amplitude A, this is on the oscillating 
side close to the bifurcation, \Cg — C c \/C c ~ (5/ A) 2 . Two 
modes of signal detection are possible: (i) For transient 
stimuli short compared to r the system operates at Cs- 
The amplitude |xi| shows the characteristic nonlinear re- 
sponse discussed above, (ii) For stimuli sustained over 
longer times, self-tuning maintains \x±\ ~ S constant for 
different stimulus amplitudes. This effect of the feedback 
represents a perfect adaptation mechanism. However, in 
the presence of noise, phase-locking of the response (to 
be discussed later) occurs as soon as an external stimulus 
is applied, and this can be detected. 



II. MODEL 

Mechanosensitivity and self-tuning mechanism 
provided by transduction channels. We demon- 
strate the general principles introduced above by devising 
a specific model for the amplification of acoustic stim- 
uli by hair bundles in non-mammalian vertebrates. A 



1 Self-tuning to a coexistence point has been discussed pre- 
viously in certain dynamic first-order transitions 



schematic representation of a hair bundle is shown in Fig. 
2. It consists of several stereocilia and a single kinocil- 
ium |l^,^0|. Transduction of hair bundle deflection to 
a chemical signal occurs via channels located near the 
tip of each stereocilium. Tip links which connect neigh- 
boring stereocilia are believed to be the gating springs of 
the transduction channels. If the hair bundle is deflected, 
tension in the tip links triggers the opening of the chan- 
nels. The subsequent influx of ions (principally K + , but 
also Ca 2+ ) causes a corresponding change of the mem- 
brane potential of the hair cell which, in turn, generates 
a nervous signal. 

The mechanosensor can be used for self-tuning, as well 
as signal detection. Many physiological processes are reg- 
ulated by ionic concentrations, so it is natural to identify 
the Ca 2+ concentration with the control parameter C. 
We assume that C decreases if the transduction chan- 
nels are closed, owing to the action of pumps in the cell 
membrane. When the hair bundle is deflected by x, the 
transduction channels open with probability P (x). We 
therefore characterize the mechanosensor by the equation 



dC_ 

dt 



C 



J P {x) 



(6) 



where J is the Ca 2+ flux through open transduction 
channels. Note that this equation provides self-tuning. 
In this case it replaces the more simple but less realistic 
Eq. (^). For our numerical examples, we use a two-state 
model for the channels with 



Po(x) 



1 



1 + Ae-' x / & 



(7) 



where (I + A)^ 1 <C 1 is the probability that a chan- 
nel is open when the hair cell is quiescent and 5 is the 
characteristic amplitude of motion to which the system 
is sensitive. For a sufficiently long relaxation time r, the 
slow variation Cq of C ~ Cb(t) + C\e lu}t can be separated 
from the small-amplitude oscillations, giving 



d t C ^-— + JoPo(\xi\ 2 ) 
T 



(8) 



where P a = J Q dt P (xie luJt + X-ie~ ZUJt ) is the averaged 
probability of channel opening in the presence of oscilla- 
tions, which increases monotonically with their ampli- 
tude X\. For physically relevant parameter values, the 
system reaches a steady state close to the bifurcation 
point, independent of the initial value of C. 

Oscillations generated by molecular motors. We 
still have to specify the nature and location of the oscil- 
lator within the hair bundle. It has been suggested that 
the transduction channels might be the source of the in- 
stability |^ 12 1; or, that myosin motors in the stereocilia 
might generate the force necessary to move the bundle 
pT| , p2"| . We propose a third possibility: the kinocilium 
could vibrate using its internal dynein motors. Recently, 
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a simple physical mechanism has been proposed which al- 
lows motor proteins operating in collections to generate 
spontaneous oscillations by traversing a Hopf bifurcation 
p5|p3| . Typically, motors move along cytoskeletal fila- 
ments and elastic elements oppose this motion. In this 
case, two possibilities exist: The system either reaches 
a stable balance between opposing forces, or it oscillates 
around the balanced state. Three time scales characterize 
this behavior: The relaxation time X/K of passive relax- 
ation, where A is the total friction and K is the elastic 
modulus; and the times w-f 1 and ui^ , where u\ is the 
kinetic rate at which a motor detaches from a filament, 
and u>2 is the attachment rate. An explicit solution of 
a simple model is derived in Appendix B. We find that 
for an appropriate value of a control parameter C, which 
is related to the ratio ojx/uj^i & Hopf bifurcation occurs 
with critical frequency 



Ka 
~X 



1/2 



(9) 



which is the geometric mean of the passive relaxation rate 
K/X and the typical ATP hydrolysis rate a = u>\ + u>2- 
The above identification of C with the Ca 2+ concentra- 
tion is consistent with the fact that Ca 2+ regulates motor 
protein activity p4j . The system oscillates if the elastic 
modulus does not exceed a maximal value Am ax ~ koN, 
where ko is the crossbridge elasticity of a motor and N 
is the total number of motors. The maximal frequency, 
obtained when K = K max , can be significantly higher 
than the ATP hydrolysis rate a. 
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FIG. 2. A hair bundle consists of a single kinocilium and 
tens to hundreds of stereocilia. The kinocilium contains 
dynein motors (red). The stereocilia contain transduction 
channels (blue) which are gated by the tension in the tip links 
(cyan); the steady-state tension is maintained constant by 
adaptation motors (magenta). In our model, the kinocilium 
is the active part of the mechanoreceptor and the stereocilia 
act as the detection apparatus. Feedback control of the active 
amplificatory process involves the influx of Ca 2+ ions through 
the transduction channels; the ions diffuse into the kinocilium 
and regulate the generation of force by the motor proteins. 



Characteristic frequency of a vibrating kinocil- 
ium. The kinocilium is a true cilium containing a cylin- 
drical arrangement of microtubule doublets and dynein 
motors. Since cilia have the well-established tendency 
to beat and vibrate with frequencies from tens of Hz 
up to at least 1 kHz |^5|,^6| , the kinocilium is a natu- 
ral candidate to be responsible for the Hopf bifurcation. 
A simple two-dimensional model can be used to discuss 
the main physical properties of a vibrating cilium near a 
Hopf bifurcation p7| . In this model, motors induce the 
bending of a pair of elastic filaments separated by a dis- 
tance a (corresponding to the distance between neighbor- 
ing microtubule doublets in the axoneme). An isolated 
kinocilium of length L and bending rigidity k, fixed at 
the basal end but free at its tip, will vibrate in a wave-like 
mode with wavelength A ~ 4L. A typical displacement 
z of the motors along the filaments leads to bending of 
the filament pair and a deflection of the tip by a dis- 
tance x ~ zL/a. The elastic bending energy is of order 
(k/L 3 )x 2 ~ (k/ 'a 2 L)z 2 . Therefore, the total elastic mod- 
ulus experienced by the motors is given by 



K 



La 2 



(10) 



The viscous energy dissipation per unit time due to 
motion of the kinocilium is of order r\L(dix) 2 ~ 
(r]L 3 /a 2 )(dtz) 2 . Therefore, the total friction experienced 
by the motors can be written as 



A 



V L 3 



XopL 



(11) 



where Ao describes the dissipation within the kinocilium 
per motor and p = N/L denotes the number of dynein 
motors per unit length along the axoneme. If internal 
friction can be neglected, i.e. if L Lq = (Xqcl 2 p/rj) 1 / 2 , 
the frequency of a vibrating cilium at the bifurcation 
point is given by 



Ka \ 1 



V 



L 2 



(12) 



Using typical values Ao — 10 -9 kg/s, a ~ 20 nm, p ~ 
5 10 8 m _1 and r/ ~ 10 _3 kg/ms, we find Lq ~ 200 nm, 
shorter than typical kinocilia. Using a ~ 10 3 s _1 , ko ~ 
10 _3 N/m and n ~ 4 10~ 22 Nm 2 (the bending rigidity 
of 20 microtubules), the frequency range between 100 Hz 
and 10 kHz can naturally be spanned by changing the 
length of the kinocilium between 1 (im and 10 pm. 

The above argument neglects the contribution of the 
stereocilia to elasticity and motion. The elastic response 
of stereocilia to hair bundle deflections has been mea- 
sured |6|]2^| . It can be well described by an angular spring 
at the base of each stereocilium which contributes an 
elastic energy per stereocilium of the order of k s (x/L) 2 , 
where k s is an angular elastic modulus. The kinocilium 
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length L in a hair bundle is approximately inversely pro- 
portional to the number N s of stereocilia Jl9|], and we 
write N s ~ where l s ~ 10~ 4 m is the total length of 
stereocilia. With this assumption, we find an additional 
elastic modulus 



contributed by the stereocilia to K. The contribution of 
stereocilia to the friction coefficient A can be estimated 
as 

\ s ~ V l s L 2 /a 2 . (14) 

The measured value of k s [p8| indicates that K s domi- 
nates the contribution to K given byEq. ©. Similarly, 
since l s > L, X s should dominate friction. In this case, 
we expect 

*"(ssr <i5) 

and the range of frequencies is somewhat reduced. 

If the stiffness of the ensemble of stereocilia greatly 
exceeds that of the kinocilium, a new situation arises. 
Since the kinocilium is attached to the stereocilia at its 
tip ]2(3] , movement of the tip is strongly reduced and the 
kinocilium preferentially vibrates in a mode with wave- 
length A ~ 2L for which the relation given by Eq. jl^ ) 
is again valid. 

Adaptation motors. Finally, we note that in order 
to obtain a robust self-tuning, the feedback mechanism 
must be sensitive only to the oscillation amplitude x\ 
and not to a stationary displacement xq', we assumed in 
Eq. (|J) that P a is independent of Xq. The transduction 
mechanism of stereocilia and their tension-gated chan- 
nels does indeed have this property: It is well known 
that an ATP-dependent adaptation mechanism p3] ex- 
ists which removes the dependence of the channel cur- 
rent on a constant displacement xq . It is widely believed 
that this adaptation involves the motion of myosin mo- 
tors, which maintains constant the steady-state tension 
in the tip links that control the transduction channels 
. Therefore, the stereociliar transduction mechanism 
has precisely the required properties to be used as a feed- 
back signal for self-tuning the bifurcation. 

III. SIMULATION 

The above analysis does not include the effects of noise. 
Brownian motion is one source of fluctuations in the 
movement of the hair bundle pl[ . Another source of 
noise is caused by stochastic fluctuations in motor pro- 
tein force as dynein molecules bind to, and detach from 
the microtubules in the kinocilium. In order to investi- 
gate the consequences of this randomness, we performed 



a Monte Carlo simulation of the two-state model de- 
scribed in |l5| ], using a realistic number of motor pro- 
teins. The motors, when attached, experienced a po- 
tential of amplitude U and period I (the corresponding 
crossbridge elasticity is fco = U/l 2 ). The attachment rate 
was constant and independent of position, u>2 = ol. De- 
tachment was localized to a region of width 0.1 1, centred 
on the potential minimum and the detachment rate u>\ 
was regulated by the Ca 2+ concentration C. We chose 
k>i/u>2 = (C q /C) 3 , where C q is the steady-state Ca 2+ 
concentration in a quiescent hair cell; the precise func- 
tional dependence is unimportant as long as oj\/u>2 de- 
creases monotonically with increasing C in a fairly sen- 
sitive way. We simulated systems with N = 1000-4000 
motors, representing hair cells with different kinocilium 
length and stereociliary number. 



3 r 




1000 2000 

at 



FIG. 3. Self-tuning of a hair bundle. If the Ca concen- 
tration C within the cell is artificially high, the hair bundle 
is initially quiescent; but as Ca 2+ ions are pumped from the 
cell, it gradually begins to oscillate with small amplitude. If C 
is suddenly artificially lowered, the hair bundle becomes un- 
stable and executes high-amplitude spontaneous oscillations. 
Since these movements open the transduction channels, the 
influx J of Ca 2+ ions increases; the consequent change in C 
regulates the motors and diminishes the amplitude of the os- 
cillation until it almost disappears. In the steady state, the 
bundle executes self-tuned critical oscillations (STCO). For 
this simulation, we assumed that the probability of transduc- 
tion channels opening was Po{z) = (1 + We- 20z/l )-\ where 
z is the motor displacement, and that the time constant for 
equilibration of the Ca 2+ concentration was r = 1000/a. 

Self-tuning and the characteristic frequency of 
spontaneous oscillations. The self-tuning of a hair cell 
to the vicinity of the bifurcation, where small-amplitude 
spontaneous oscillations occur, is demonstrated in Fig. 
3. When a change of the internal Ca 2+ concentration 
is imposed, the system is transiently perturbed; but af- 
ter an interval of time of order r, it returns to the same 
steady state. Hair cells with different numbers of motors 
acquire different internal concentrations of Ca 2+ , in or- 
der to adjust the motor detachment rate in such a way 
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that the system approaches the critical point. The spon- 
taneous oscillations of three different hair bundles are 
shown in Fig. 4. Note that the characteristic frequency 
is approximately proportional to the inverse-square of the 
number of motors, as expected from Eq. (|l2|), and that 
it can exceed the typical ATP cycle rate a when the total 
number of motors is small (short kinocilium). All three 
hair cells execute spontaneous oscillations with a similar 
amplitude, as expected from our arguments above. How- 
ever, the noise introduces a significant new effect: The 
oscillations are irregular. The incoherence of the phase 
of the oscillation is evident in the Fourier transform of 
the displacement, which exhibits a broad peak centred 
on the characteristic frequency. 
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FIG. 4. Self-tuned critical oscillations of systems compris- 
ing N = 1000, 2000 and 4000 motors. In each case, the os- 
cillation has small amplitude z ~ 0.1 1 and is irregular, as 
can be seen from the broadness of the Fourier spectrum. For 
these simulations we chose parameters ko/Xa — 4 W 7 N~ 3 
and K/Xa = 4 10 13 iV -4 , which correspond with the scaling 
dependence of the spring constant and friction coefficient on 
kinocilium length in Eqs. dicf-flll) , and also with the order of 
magnitude estimates of the physical parameters, given in the 
text. 

Dynamic response to a tone at the characteris- 
tic frequency. The response of a self-tuned hair bun- 
dle to a sinusoidal force with a frequency approximately 
equal to the bundle's characteristic frequency is illus- 
trated in Fig. 5. For weak stimuli, the amplitude of the 
oscillation does not increase with the amplitude |/i| of 
the applied force; this is because the small response to 
the stimulus is masked by the noisy, spontaneous motion. 
Instead, the phase of the hair-bundle oscillation becomes 
more regular; as it does so, a peak emerges from the 
Fourier spectrum at the driving frequency. The height of 
the peak grows approximately as |/i| 1,/3 for intermediate 
values of |/i|, and approximately linearly for very weak 
stimuli, and also for very strong stimuli (for which the 



response is essentially passive). Thus the Fourier compo- 
nent of the hair-bundle displacement at the driving fre- 
quency responds to the stimulus in the generic manner 
discussed above. 
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FIG. 5. Response of a system with TV = 2000 motors to a si- 
nusoidal force at frequency v — 0.3a, close to the hair bundle's 
characteristic frequency. Curves are marked by the dimen- 
sionless amplitude of the force, \fi\/fmot, where f mo t = U/l 
is the force produced by a motor molecule. Note that a force 
equal to that of a single motor is sufficient to elicit a response 
in the Fourier spectrum. Curves shaded red are the responses 
of an equivalent, passive hair bundle (i.e. a bundle with iden- 
tical mechanical properties but no force generators). 

Phase-locking of the nervous signal. The flow of 
ions through the transduction channels depolarizes the 
cell membrane which, in turn, opens voltage-gated chan- 
nels at the base of the hair cell and generates a synaptic 
current Jl[ . This sequence of events happens fast enough 
for variations at the synapse faithfully to reflect the hair- 
bundle motion at frequencies below 1 kHz. Information 
about the auditory stimulus is subsequently passed along 
the auditory nerve in the form of a spike train. Simpli- 
fying this transduction process, we assume that a spike 
is elicited whenever the transduction channel current J 
passes a threshold value of 0.5 J . The resulting nervous 
response is shown in Fig. 6. In the absence of a stimu- 
lus, the self-tuned critical oscillations of the bundle cause 
the nerve to fire stochastically, at a low rate. When 
a weak sinusoidal stimulus is applied at the character- 
istic frequency, the firing rate does not increase above 
the spontaneous rate, but the spike train becomes de- 
tectably phase- locked to the stimulus. The degree of 
phase-locking increases rapidly as the amplitude of the 
stimulus increases, reflecting the growing regularity of 
the hair-bundle motion. It is only when the neural re- 
sponse is almost completely phase-locked that the firing 
rate begins to rise. Eventually, for strong stimuli, the 
spike rate saturates at the stimulus frequency. This be- 
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havior is strikingly similar to that which is observed ex- 
perimentally. In particular, it is well known that the 
threshold for phase-locking in the auditory nerve fiber 
is 10-20 dB lower than the threshold at which the firing 
rate begins to rise [p2 33 1. 
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FIG. 6. (a) Histogram of the phase <f> of the driving force 
at the instants when nervous spikes are generated, for var- 
ious values of the force amplitude (marked by the value of 
\fi\l fmot)- (b) Degree of phase- locking (blue; expressed as 
the vector strength J d<j> p(^)e"*) and firing rate (red), as a 
function of the driving force. The frequency selectivity of the 
bundle can be appreciated by comparing the neural response 
to a stimulus with frequency close to the characteristic fre- 
quency, v = 0.3a (•), with the response to other frequencies, 
v = 0.15a (v) and v = 0.6a (A). 



IV. DISCUSSION 

The self-tuned critical oscillator which we have intro- 
duced as a system for signal detection has characteris- 
tic amplification properties which are generic and do not 
depend on the choice of a specific model. The oscilla- 
tion frequency, however, depends on the physical mech- 
anism involved. In our model of non-mammalian hair 
cells, the frequency is determined by the geometry of the 
hair bundle. A simple morphological gradient along the 
basilar membrane would endow the ear with the abil- 
ity to analyze a wide range of frequencies. Experimen- 
tally, it is well established that the height of hair bun- 
dles progressively increases along the cochlea, and that 
concurrently the characteristic frequency of the hair cells 
declines [ p^p^[ |. Our proposition that the kinocilium is 
likely to play an active role in non-mammalian vertebrate 
hair cells suggests experiments which study the motility 
of the kinocilium and its potential for generating oscilla- 
tions. If the kinocilium is the only source of oscillations, 
its removal should suppress hair-bundle vibrations. If 



the self-tuning mechanism is removed instead (by cut- 
ting the tip links, for example) our model predicts that 
the kinocilium should exhibit stronger spontaneous oscil- 
lations. Careful control of the extracellular ion concen- 
trations, such as Ca 2+ would be essential in such exper- 
iments. 

The simulations of our model reveal how a hair bundle 
can achieve its remarkable sensitivity to weak stimuli. By 
profiting from the periodicity of a sinusoidal input, and 
measuring phase-locking rather than the amplitude of re- 
sponse, the mechanosensor can detect forces considerably 
weaker than those exerted by a single molecular motor (if 
the bundle were a simple, passive structure, its response 
to such forces would be smaller than its Brownian mo- 
tion). An important implication of this detection mech- 
anism is that even though the hair cell selects a certain 
frequency, the signal must still be encoded by the interval 
between spikes elicited in the auditory nerve. Paradoxi- 
cally, the stochastic noise caused by the motor proteins 
serves a useful purpose. It ensures that the self-tuned 
critical oscillations of the hair bundle are incoherent, so 
that the pattern of spontaneous firing in the nerve is ir- 
regular. Against this background, the regular response 
to a periodic stimulus can easily be detected.^ 

Another beneficial feature of noise arises from the fact 
that weak stimuli do not increase the amplitude of oscil- 
lation above the spontaneous amplitude. Thus the Ca 2+ 
concentration remains constant, the hair bundle stays in 
the critical regime, and active amplification can be sus- 
tained indefinitely. Stronger stimuli cause the system to 
drift away from the critical point, so that the degree of 
amplification diminishes over time. It is well known that 
both the perception of loudness |?5| and the firing rate 
of the auditory nerve fl36| decrease over a period of a few 
seconds, when a stimulus of moderate intensity is main- 
tained. This phenomenon, which is usually referred to 
as 'adaptation', is consistent with our self-tuning mech- 
anism. Following the sustained presentation of a loud 
stimulus, the spontaneous firing rate diminishes and the 
threshold to weak stimuli is augmented |3(|. Such 'fa- 
tigue' is also accounted for by our self-tuning mechanism. 

Since self-tuning positions the system slightly on the 
oscillating side of the critical point, self-tuned criticality 
provides a natural explanation for otoacoustic emissions 
In its normal working state, the inner ear would 
generate faint sounds with a broad range of frequencies. 
If the feedback mechanism were to fail in certain cells, 
the spontaneous oscillations could become large enough 
for distinct tones to be emitted. 

A self-tuned Hopf bifurcation is ideal for sound detec- 
tion because it provides sharp frequency selectivity and a 



2 The benefits of noise have been discussed in a variet; 
situations such as those involving stochastic resonance 
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nonlinear gain which compresses a wide range of stimulus 
intensities into a narrow range of response. We therefore 
believe that the concept applies to all vertebrate hear- 
ing systems, and potentially to other mechanoreceptor 
systems. Kinocilia are typically absent in mammalian 
cochlea, and we suggest that their force-generating role 
has been assumed by the outer hair cells. Self-tuning of 
these motile cells might be realized by a mechanism sim- 
ilar to that presented here, using transduction channels 
in their hair bundles. It could, however, work very dif- 
ferently; for example, it might involve feedback from the 
inner hair cells, via the efferent nervous system. 

Tuning to the proximity of a critical point is likely to 
be a general strategy adopted by sensory systems. Simple 
molecular receptors Q , as well as the physiological sen- 
sors of higher organisms, can enhance their response to 
weak stimuli in this way. We propose that the physics of 
self-tuned criticality is the 'central science of transducer 
physiology' spoken of by Delbriick |3j| . 



change the basic properties. In this situation, the rela- 
tion between x and / can be expressed as a systematic 
expansion of the force amplitudes f n in the amplitudes 
x n . 

fk = ^ Xi + F^^XlXm 

^klmn XlXmXn 

+ 0(x i ) , (A3) 

where the expansion coefficients Fjfjj k are symmet- 
ric with respect to permutations of the indices k\..k n . 
The limit cycle solutions are invariant with respect to 
translations in time t — > t + At. Under these transfor- 
mations the amplitudes change e m 4 and 
In -> f n e inuAt . Inspection of Eq. (|a|) shows that the 
time translation symmetry allows only for those terms 
•^fc'fei k X kx--Xk n for which k = k\ + ... + k n . For all 
other cases F kki k must vanish which significantly re- 
stricts the number of terms. 
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APPENDIX A: GENERIC BEHAVIOR AT A 
HOPF-BIFURCATION 

1. Nonlinear relation between periodic stimulus and 
displacements 

We are interested in the response x(t) of a nonlinear 
system to a periodic stimulus force fit). If only one fre- 
quency v = lu/2it is present we use the Fourier expansions 



/(f) = £ / n e* 

71 — — OG 
OO 

x(t) = 2J x ne l 



(Al) 
(A2) 



n— — oo 



where the complex coefficients x n and /„ obey x n = x*_ n 
and /„ = fl n . This representation implies that we focus 
on the limit cycle solution and ignore all transient relax- 
ation phenomena. We consider the class of systems for 
which the force at a given time depends in a nonlinear 
way on the history of the displacements x(t) alone; as 
we will discuss in section D more complex cases do not 



2. Hopf bifurcation 

The nonlinear system exhibits spontaneous oscillations 
and a Hopf-bifurcation if nontrivial solutions to Eq. (A3) 
with x n ^ exist in the case where all fk = 0, i.e. if no 
stimulus force is applied. Without loss of generality, we 
consider here an instability of the mode x\. In this case, 
the dominant terms allowed by symmetry read (fk = 0) 







2-^1: 2 _ ± X-iX2 



F^X! 



6 - p lfl,l -1^1^-1 



F^x 2 



'^F2i\x 1 



6-Pll 2 -2 X 2 X ~2 X 'L 



Eq. (|A5|) determines x 2 ^ — 2(i ? 21 i/-^22 ) x i- Inserting 



(A4) 
(A5) 



this relation in Eq. (A4), we obtain to lowest order 



~ Axi +B\x 1 \ 2 x 1 



where A = F$ and B : 



1,1,1, 



7(2), 
211 J 



p(2) 
1,2.- 



(A6) 

l/ r 22 ■ 



The coefficients A(ui, C) and B(ui, C) are complex and 
in general depend on frequency to and a control parameter 
which we denote by C. A Hopf bifurcation occurs at a 
critical point C = C c at which A vanishes for a frequency 
u) c , i.e. A(lo c ,C c ) = 0. This can be demonstrated as 
follows: A spontaneously oscillating solution satisfies 



A 
B 



(A7) 



Note, that such a solution can only exist if A/B is real 
and negative. At the bifurcation point, .4 = and A/B 
is therefore real for lj = u> c , however the corresponding 
amplitude \xi\ 2 vanishes. In the vicinity of this point we 



expect to find solutions with finite amplitude, 
the expansion 



We use 
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A{w, C)~(u- lo c )A 1 + (C - C C )A 2 



(A8) 



where A\ and A 2 are complex coefficients and we neglect 
higher order terms. Spontaneous oscillating solutions ex- 
ist only if A/B is real. This condition is satisfied for a 
particular frequency uj = uj s with 



UJ S 



Im{A 2 /B) 
Im{A x /B) 



{C c - C) 



(A9) 



The ratio —A/B at this frequency uj s changes sign for 
C = C c ; here we assume without loss of generality that it 
is positive for C < C c . In this case, the system oscillates 
s pon taneously with an amplitude which according to Eq. 
(|A7|) behaves as |xi| 2 = A 2 (C C - C)/C c , where 



A 



C c [Re{A 2 /B) 



Re(A 1 /B) 



Im{A 2 /B) 
Im{A 1 /B) 



(A10) 



is a typical amplitude. We have thus demonstrated that 
Eq. ( |A6| ) characterizes a Hopf-bifurcation if the complex 
coefficient A vanishes at a critical point C c for a critical 
frequency uj c . 



3. Amplified response to sinusoidal stimuli 



If a sinusoidal stimulus f(t) — j\e lu)t 
which all /„ with n ^ ±1 vanish, Eq. (|A 

h ~Aei +B\ Xl \ 2 Xl . 



|_ f_ ie -i^t ) for 
j) becomes 

(All) 



We consider a system that is tuned exactly to the bifurca- 
tion, C = C c . In this situation spontaneous oscillations 
do not occur and A — (uj — u> c )Ai. If the imposed fre- 
quency is equal to the critical frequency uj = uj c , the 
coefficient A vanishes and we can solve Eq. (All) for 
|a;i| to find the nonlinear response 



\B\ 



l/3| A |l/3 



(A12) 



as a function of the force amplitude |/i|. This behavior 
represents an amplified response with a gain 



Fil 
l/il 



l/i 



-2/3 



(A13) 



that becomes arbitrarily large for small forces. If the 
frequency uj is different from uj c , this nonlinear response 
still holds as long as the linear term in Eq. (All) is small 
compared to the cubic term and can be neglected. This 
is the case if |£i| 2 > \A/B\ = \uj - u c \\Ai/B\. There- 
fore, the nonlinear regime characterized by Eq. (A12) 
holds for sufficiently large force amplitudes, |/i| 3> \(uj — 
uj c )Ai \ z / 2 /\B\ X / 2 , or if the frequency is sufficiently close 
to the critical frequency, \uj - uj c \ < |/i| 2/3 |#| 1/3 /l^i|- 

If the frequency mismatch \u> — lu c \ becomes large, or 
if forces are small, a new regime occurs for which 



the linear term in (All) dominates. In this regime, the 
response is linear, 



l/i 



\(uj - uj c )Ai\ 



(A14) 



and the gain is constant. This is a passive response if the 
stimulus frequency is too far from the critical frequency. 



4. Additional remarks 

The above derivation is based on an expansion ( |A3| ) in 
the displacements x n . This excludes some nonlinearities 
in the force which can lead to additional nonlinear terms 



in Eq. (All). The most general form of Eq. (All) is 



fi ~ Ax i 



B\x 1 \ 2 x 1 
£W\ 2 h- 



- Cxi\h\ 



+ Vx- 
<?l/i| 2 /: 



if? 



(A15) 



However, for small forces /i and small amplitudes x±, 
the results derived above are not affected. The regime 
of nonlinear reponse |/i| ~ | a^i | 3 , as well as the linear 
response regime |/i| ~ |xi| still exist. If |/i| ~ |xi|, the 
nonlinear terms in fi renormalize the third order term, 
which in this regime is negligable. If |/i| ~ |a^i| 3 , the 
nonlinear terms in fi are of even higher order and can 
be neglected. 



APPENDIX B: OSCILLATIONS GENERATED BY 
MOLECULAR MOTORS 

1. Two state model 

The two state model describes force-generation as a re- 
sult of transitions between two states, a bound state and 
a detached state of a motor and its track filament. The 
interaction between a motor at position z along the fila- 
ment in states 1 and 2 is characterized by two periodic 
potentials Wi(z) = Wi(z + I) and W 2 {z) = W 2 (z + I) 
where / is the period. We introduce the relative position 
£ = z mod I with respect to the potential period. De- 
tachment and attachment rates are denoted k>i(£) and 
w 2(£)j respectively. Oscillations can occur in this model 
if a large number N of motors move collectively against 
an external elastic element of modulus K . 

We introduce the probability Pi (£ ) and P 2 (£) of finding 
a motor bound at position £ in state 1 or 2, which satisfy 
the normalization condition 



/ d£(Pi+P 2 ) 
Jo 



= 1 



(Bl) 



For a large number of motors collectively moving with 
the same velocity v the dynamic equations read 
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d t Pi + vd^Pi = -ujiPi + w 2 P 2 
d t P% + vd^P 2 = u)\P\ - UJ2P2 



(B2) 
(B3) 



The velocity v is determined by the force-balance condi- 
tion 

f = Xv + Kz + N [ d£(Pidt:Wi + P 2 d i W 2 ) (B4) 
J 

where A is a friction coefficient describing the total fric- 
tion and z is the displacement of the motors, dtz = v. For 
an incommensurate arrangement of motors with respect 
to the track filament and a large number N of motors, 
PliO + = 1/' and the equations of motions sim- 

plify: 

d t P + vd i P=-(LU 1 +UJ 2 )P + LU 2 /l , (B5) 

where we denote for simplicity P(£) = Pi(0- 

We discuss a simple choice for the potentials and tran- 
sition rates for which the Hopf bifurcation is easy to de- 
termine analytically. We consider the potential 



Wx(£) = Z7cos(2tt£/0 



(B6) 



with amplitude U , and the potential W 2 to be constant. 
The transition rates are chosen to be periodic functions 



wi(0 =/?-/?cos(2</0 
wa(0 =a-0 + 0cos(27rf/O 



(B7) 
(B8) 



parameterized by two coefficients a and j3. With this 
choice, 



+ uj 2 (0 = a 



(B9) 



is constant and the fact that u>\ and lo 2 are positive re- 
stricts (3 to the interval < (3 < a/2. 



2. Linear response function 

In order to determine the linear coefficient A which 
determines the stability of the system, we look for small 
amplitude oscillations close to the resting state with v = 
0. We write 



P~ Pa + Pl e wt 
f - fie luJt 



z ~ Z\e 



(BIO) 
(Bll) 
(B12) 



wher e pp = uj 2 /al. To linear order in z\, we find from 
Eq. © 



Pi 



IUJZ\ 



d x po 



iuj + a 

The corresponding force is given by 



(B13) 



fi ^ Az x 



(B14) 



with 



A = iwA + K + x , (B15) 
where the active response x of the motors is given by 



X =-N ^-—d^Wx (B16) 
J iuj + a 



For the choice of Eq. ( p36| ) and (BS) the integral can be 
calculated and we obtain 

A(C, u>) = iu\ + K- Nk C iUJ { a + { "{f . (B17) 

Here, we have introduced the dimensionless control pa- 
rameter C = 2ir 2 f3/a with < C < it 2 and the cross- 
bridge elasticity fco = U/l 2 of the motors. 



3. Hopf bifurcation 

A Hopf bifurcation occurs if there is a pair of values 



(C, ui) for which A as given by Eq. (B17) vanishes. Such 

(B18) 



a point indeed exists. For the critical value 
Act + K 



C c 



Nk 



the bifurcation occurs for the critical frequency 



1/2 



(B19) 



The critical frequency is bounded by the fact that C c < 
ir 2 . The maximal frequency occurs for the maximal pos- 
sible value of K 

K max = Nk Q ir 2 - \a (B20) 

for which C c = ir 2 . This frequency is given by 



1/2 



(B21) 



Note, that the maximal frequency can be significantly 
higher than the typical rate a of the chemical cycle. 
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